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Re-analysis of the solar phase curves of the icy Galilean 
satellites demonstrates that the quantitative results are depen- 
dent on the single particle scattering function incorporated into 
the photometric model; however, the qualitative properties are 
independent. The results presented here show that the general 
physical characteristics predicted by a Hapke model (B. Hapke, 
1986, Icarus 67, 264-280) incorporating a two parameter double 
Henyey-Greenstein scattering function are similar to the pre- 
dictions given by the same model incorporating a three parame- 
ter double Henyey-Greenstein scattering function as long as 
the data set being modeled has adequate coverage in phase 
angle. Conflicting results occur when the large phase angle 
coverage is inadequate. Analysis of the role of isotropic versus 
anisotropic multiple scattering shows that for surfaces as bright 
as Europa the two models predict very similar results over 
phase angles covered by the data. Differences arise only at 
those phase angles for which there are no data. The single 
particle scattering behavior between the leading and trailing 
hemispheres of Europa and Ganymede is commensurate with 
magnetospheric alterations of their surfaces. Ion bombardment 
will produce more forward scattering single scattering functions 
due to annealing of potential scattering centers within regolith 
particles (N. J. Sack et ai, 1992, Icarus 100, 534-540). Both 
leading and trailing hemispheres of Europa are consistent with 
a high porosity model and commensurate with a frost surface. 
There are no strong differences in predicted porosity between 
the two hemispheres of Callisto, both are consistent with model 
porosities midway between that deduced for Europa and the 
Moon. Surface roughness model estimates predict that surface 
roughness increases with satellite distance from Jupiter, with 
lunar surface roughness values falling midway between those 
measured for Ganymede and Callisto. There is no obvious 
variation in predicted surface roughness with hemisphere for 
any of the Galilean satellites. * 1997 Academic press 
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INTRODUCTION 

The photometric properties of a planetary surface are 
determined by the physical structure and texture of the 
surface regolith and its material. Photometric modeling 
attempts to describe and place constraints on such proper- 
ties as roughness, porosity, and regolith grain size distribu- 
tion, while also describing and placing constraints on such 
material properties as the average single scattering albedo 
and single scattering functions. All these properties are 
influenced or determined by the geological processes that 
formed and shaped the planetary surface. Therefore under- 
standing the photometrically determined physical structure 
of a planetary reolith provides supporting evidence for 
many geologic processes such as micrometeorite and ion 
bombardment. Photometric modeling can also assist in 
interpretations of spectral data, where such regolith prop- 
erties as grain size and compositional mixing can influence 
the shape of reflectance spectra. This makes photometric 
modeling a useful tool in comparative planetology. 

Background 

The solar phase curves of the icy Galilean satellites have 
been modeled using Hapke’s theory (Hapke 1981, 1984. 
1986) by several investigators (Buratti 1985, 1991, Dom- 
ingue et ai 1991, Helfenstein 1986). Buratti (1985) exam- 
ined Europa’s phase curve based on Voyager images using 
an early version of Hapke’s theory. She did not separate 
her data into leading versus trailing hemispheres, but ex- 
amined the surface of Europa as a whole. Buratti (1985) 
concluded from her analyses that the optically active por- 
tion of Europa's surface is highly compact and the regolith 
particles show nearly isotropic scattering behavior. Dom- 
ingue et ai (1991) expanded on her work by combining 
telescopic observations with the Voyager image data to 
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construct phase curves that included opposition data (that 
were not available to Buratti) and by separating the data 
into leading and trailing hemisphere subsets to examine 
any potential hemispherical dichotomies. Their modeling 
of the new opposition data with an updated version of 
Hapkie's theory (Hapke 1986) predicts that the optically 
active portions of the regolith on both hemispheres of 
Europa are very porous, consistent with a frost layer. Dom- 
ingue et ai (1991) found that their results were consistent 
with the regolith particles of the leading hemisphere being 
strongly forward scattering and the regolith particles of 
the trailing hemisphere having both forward and backward 
scattering components to their scattering function. The 
differences in the particle scattering behavior found by 
Buratti (1985) and Domingue etai (1991) can be attributed 
to their choices in the single particle scattering function 
incorporated into Hapke’s model. The different compac- 
tion states found by these two studies are due to the inclu- 
sion of the telescopic data, which measures the opposition 
effect in more detail and to a smaller phase angle than the 
Voyager data alone. Helfenstein (1986) applied Hapke’s 
model to individual geologic terrains on Ganymeade using 
Voyager image data. The terrains on Ganymede can be 
divided into two principle units, bright and dark. Hel- 
fenstein (1986) found that on average the photometric 
properties, as deduced from Hapke’s model, for the dark 
terrains show a greater macroscopic roughness, a higher 
state of compaction, a more uniform single scattering be- 
havior, and a lower single scattering albedo than the bright 
terrains. Buratti (1991) examined the disk-integrated phase 
curves of both Ganymede and Callisto using a combination 
of telescopic observations and measurements from the 
Voyager images. She divided her data sets into leading and 
trailing hemispheres in order to examine any hemispherical 
dichotomies in the photometric surface characteristics of 
these two satellites. Buratti (1991) found that Ganymede 
displays no hemispherical differences in surface roughness, 
compaction state of the optically active layer, or scattering 
behavior of the regolith particles. She attributed these re- 
sults to the balance of magnetospheric interactions on the 
trailing hemisphere with micrometeoritic erosion on the 
leading hemisphere. Buratti (1991) did find hemispherical 
differences on Callisto where the optically active layer of 
the trailing hemisphere’s regolith is more compact than the 
leading hemisphere’s regolith and the leading hemisphere 
regolith particles are more backscattering than those on 
the trailing hemisphere. 

The research presented here is a reexaminationof the 
disk-integrated phase curves of all three of the icy Galilean 
satellites. The phase curves examined are constructed from 
combinations of telescopic observations and Voyager im- 
age measurements, following the methods of Buratti (1991) 
and Domingue et ai (1991). New, previously unpublished, 
telescopic observations by Lockwood and Thompson (pri- 


vate communication) are included for Callisto. All of the 
studies discussed above assumed the multiply scattered 
component of surface reflected light is isotropic. This is not 
a valid assumtion for very bright objects, such as Europa 
(Mishchenko 1994, Goguen 1996, Verbiscer and Hel- 
fenstein 1997). Thus the purpose of this study is twofold: 
first, it examines the effect of the choice of which single 
particle scattering function is incorporated into the basic 
Hapke model where isotropic multiple scattering is as- 
sumed. The goal is to examine if the hemispherical differ- 
ences in surface structure and texture measured by Buratti 
(1991) and Domingue et ai (1991) are real or an artifact 
of the choice in single particle scattering function. Second, 
this paper explores the differences between Hapke model- 
ing results for isotropic and anisotropic modeling of the 
multiple scattering component for bright objects such as 
Europa and Ganymede. 

Model 

An integral part of any photometric theory which at- 
tempts to model scattering from a planetary regolith is the 
incorporation of a single particle scattering function. The 
single particle scattering function’s role is to describe how 
an idealized, average particle or grain from a planetary 
regolith scatters light. Such knowledge is an important 
verification tool for various geologic processes. In Hapke’s 
model it is assumed that all regolith particles or grains 
are large compared to the wavelength of incident light. 
McGuire and Hapke (1995) studied the changes in the 
single particle scattering function as a function of particle 
shape, texture, and clarity. Using 1-cm (in diamter) sized 
artificial resin particles McGuire and Hapke measured the 
scattering functions as shape, particle surface texture, and 
density of internal scattering material were systematically 
varied for incident wavelengths of 448, 554, and 690 nm. 
They found that to model the resulting scattering phase 
curves of their artificial particles a two parameter single 
particle scattering function was needed and that a two 
parameter double Henyey-Greenstein function modeled 
their results better than a two parameter Legendre polyno- 
mial. Hartman and Domingue (1997) compared the quality 
of fit between a two parameter double Henyey-Greenstein 
function versus a three parameter double Henyey- 
Greenstein function to the McGuire and Hapke (1995) 
scattering data. The goal of their research was to establish 
which type of scattering function was essential for ob- 
taining physically meaningful results from Hapke’s model. 
Helfenstein et ai (1997) also fit the McGuire and Hapke 
(1995) scattering data with a three parameter Henyey- 
Greenstein function to examine the relationship between 
scattering albedo and single scattering function. The dou- 
ble Henyey-Greenstein functions are defined and dis- 
cussed under Analytical Methods. In the course of their 
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study Hartman and Domingue (1997) found that the choice 
of single particle scattering function influenced the values 
of the non-opposition Hapke parameters (the opposition 
Hapke parameters are those which describe the opposition 
effect, all others are non-opposition parameters). 

In the first set of analyses presented here the phase 
curve data are modeled twice, once incorporating a two 
parameter double Henyey-Greenstein function and then 
a second time incorporating a three parameter double Hen- 
yey-Greenstein function into Hapke’s model. This differs 
from the work done by Buratti (1985, 1991) in that she 
used a single parameter Henyey-Greenstein function. 
Domingue et ai (1991) used a two parameter double Hen- 
yey-Greenstein function; therefore, this work will expand 
on their findings with the three parameter scattering func- 
tion application. The goal of this initial study is to examine 
if the hemispherical dichotomies observed in the phase 
curve analysis of these satellites is a function of the choice 
in single particle scattering function. This first set of analy- 
ses assumes that the multiply scattered component is iso- 
tropic. 

The second set of analyses investigates the differences 
in modeling results if anistropic multiple scattering is incor- 
porated into the model. For the second part of this study 
the phase curves were reanalyzed using a three parameter 
double Henyey-Greenstein function where it is assumed 
that the anisotropy in the singly and multiply scattered 
components can be described by the same scattering func- 
tion (Verbiscer 1991). The anisotropic multiple scattering 
analyses use the method described by Verbiscer (1991). 

DATA 

Earth- Based Observations 

The Earth-based telescopic observations used in this 
study are taken from Domingue et ai (1991) for Europa 
(which are the observations of G. W. Lockwood and 
D. T. Thompson), from Millis and Thompson (1975) for 
Ganymede, and from Lockwood and Thompson (private 
communication) for Callisto. The Europa and Callisto tele- 
scopic observations were taken in the Stromgren y and b 
filters, which have effective wavelengths of 0.55 and 0.47 
/xm, respectively. The Ganymede observations from Millis 
and Thompson (1975) were taken in the V and B filters, 
which correspond to 0.55 and 0.44 /xm. respectively. All 
telescopic data were corrected for rotational phase varia- 
tions. The rotational phase angle ranges for the leading 
and trailing hemispheres of all three satellites are listed in 
Table I. The broad rotational phase angle ranges make 
any measured leading/trailing hemisphere dichotomies de- 
pendent on the rotational phase curve corrections applied 
to the data sets. The leading hemisphere of Europa was 
corrected to 95° rotational phase and the trailing hemi- 
sphere was corrected to 295° rotational phase using the 


TABLE I 


Satellite 

Data 

source 

0.47-/xm 

lead 

0.47-/xm 

trail 

0.55-/ mi 
lead 

0.55 -/xm 
trail 


Satellite phase 

angle coverage 


Europa 

Telescopic 

0.2- 11.3 

0.2- 11.3 

0.2- 11.3 

0.2- 11.3 


Voyager 

3.8-119.0 

3.0-119.0 

3.7-104.4 

2.9- 94.3 

Ganymede 

Telescopic 

0.6- 11.4 

1.4- 11.2 

0.4- 11.4 

1.4- 11.2 


Voyager 

4.8-120.3 

20.0-124.0 

4.9- 32.2 

20.0- 38.7 

Callisto 

Telescopic 

0.4- 11.1 

0.6- 11.2 

0.4- 11.1 

0.6- 11.2 


Voyager 

8.3- 30.4 

13.9-140.8 

8.3- 51.7 

14.0-107.4 


Satellite rotational angle coverage 


Europa 

Telescopic 

31-126 

213-302 

31-126 

213-302 


Voyager 

7-165 

207-351 

11-162 

207-346 

Ganymede 

Telescopic 

12-151 

213-337 

12-151 

213-337 


Voyager 

17-176 

185-359 

17-176 

186-360 

Callisto 

Telescopic 

32-133 

226-312 

32-133 

226-312 


Voyager 

3-174 

184-359 

3-171 

184-352 


Note. Lead corresponds lo leading hemisphere, and trail corresponds 
to trailing hemisphere. All angle ranges are in degrees. 


rotational phase curves of Domingue et ai (1991). The 
reduction of the Europa data follows that of Domingue 
et ai (1991). The leading hemispheres of Ganymede and 
Callisto were corrected to 90° rotational phase angle, 
whereas the trailing hemispheres of these two satellites 
were corrected to 270° rotational phase angle. The correc- 
tions applied to Ganymede were based on the rotational 
phase curves of Millis and Thompson (1975). The correc- 
tions applied to Callisto were based on the rotational phase 
curve observations of Lockwood and Thompson (private 
communication), Millis and Thompson (1975), and Mor- 
rison et ai (1974). Figure 1 shows the correction curve to 
255° rotational phase angle for Callisto. The rotational 
phase curves for Callisto showed the largest variation be- 
tween different data sets. The best agreement was between 
the Millis and Thompson (1975) and the Lockwood and 
Thompson Callisto data sets. For those regions of the phase 
curve where there were Lockwood and Thompson mea- 
surements, the rotational phase curve corrections based 
on their data were used. Outside of those regions the phase 
curve corrections based on the Millis and Thompson data 
were used. These corrections were used since the telescopic 
observations in the solar phase curves is predominately the 
Lockwood and Thompson measurements. The telescopic 
data, which were corrected to opposition, were also nor- 
malized to a Lambert disk of the same size at the same 
distance according to the relation 

M = A/pbmci - A/ S „n - 5 log [(s<S)/(Ar)]. 

where .v is the Earth-satellite distance (4.2 AU), S is the 
Sun-satellite distance (5.2 AU. Burns 1986), A is the Sun- 
Earth distance, and r is the satellite radius. M phinct is the 
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correction function 
Millis and Thompson 
Morrison et al. 

Lockwood and Thompson 


FIG. 1. The correction curves (solid line) used to correct the telescopic and Voyager Callisto data sets for brightness variations due to rotational 
phase. The corrections are to 255° rotational phase angle at (a) 0.47 ftm and (b) 0.55 fxm and are based on the observations of Lockwood 
and Thompson (private communication), Millis and Thompson (1975), and Morrison et al (1974), corresponding to the circles, triangles, and 
squares, respectively. 


measured magnitude of the satellite and M Sun is the magni- 
tude of the Sun at the same observing wavelength (4.378 
and 3.728 at 0.448 and 0.55 jum, respectively; Allen 1973). 

Voyager Image Data 

The Voyager images used in this study were taken in 
the clear, blue, orange, and green filters. It should be noted 


that Buratti (1991) used only clear filter images in her 
study of Ganymede and Callisto. All the images used were 
calibrated using the Planetary Image and Cartography Sys- 
tem (PICS), which incorporates the calibration corrections 
described by Danielson et al. (1991). Since the calibration 
of these images the correction factors have been updated 
(E. Eliason and G. Danielson, private communication). 
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The use of the older calibration factors produces a 7-8% 
uncertainty in the green and orange filter data and an —3% 
uncertainty in the clear and blue filter data. The calibrated 
images were converted into disk integrated values by aver- 
aging the brightness values of the pixels composing the 
satellite disk. This methods introduces a ±1% error in 
the integral data. The intensity measurements were then 
converted to magnitudes to correspond with the tele- 
scopic observations. 

The average brightness values were corrected for rota- 
tional phase angle variations using the rotational phase 
curve corrections applied to the telescopic data for the 
same satellite. The Voyager longitude system differs from 
the longitude system utilized by astronomers by 1° or 2°. 
This small difference may introduce small error bars in 
the rotational corrections. Images where the subspacecraft 
longitude fell between 0° and 180° were grouped into the 
leading hemisphere data set, and similarly, the images 
where the subspacecraft longitude fell between 180° and 
360° were grouped into the trailing hemisphere data set. 
This method poses a potential problem at high phase angles 
where the subspacecraft longitude is not centered on its 
respective hemisphere. The images used here have been 
examined for data points which have been binned in one 
hemisphere category but actually have the illuminated por- 
tion of their disk in the other hemisphere. For the images 
in this study this method bins the images correctly, where 
the illuminated portion of the disk at high phase angles is 
50% of the hemisphere it has been assigned. Table I lists the 
range of rotational and solar phase angles for the satellites. 

Images were taken by two spacecraft (Voyagers I and 
II), each having a narrow and a wide-angle camera. Within 
each spacecraft image set comparisons at similar solar 
phase angles within each filter were made between narrow' 
angle camera images and wide angle camera images. A 
correction was applied to the wide angle camera brightness 
values to match them to the narrow' angle camera bright- 
ness values if they did not coincide. The Voyager II clear 
and blue filter wide angle camera images used in the Eu- 
ropa phase curves required a correction factor of —0.13 
magnitudes, whereas the Voyager 1 1 orange and green filter 
wide angle camera images for the Europa phase curves 
required a —0.036 magnitude correction. The Ganymede 
wide angle camera images from both spacecraft required 
no corrections. The Callisto clear filter images from both 
spacecraft also required no wide angle camera corrections. 
However, the Voyager I blue filter wide angle camera 
images required a —0.198 magnitude correction, whereas 
the same Voyager II images needed a 0.145 magnitude 
correction. The corrections for the Callisto orange and 
green filter images were more complicated. The Voyager 
II wide angle camera images for these two filters required 
corrections of -0. 1 54 and — 0. 1 78 magnitudes, respectively. 
The Voyager I wide angle camera images required correc- 


tions of -0.088 and -0.178 magnitudes for these two filters, 
respectively, for the leading hemisphere images. The trail- 
ing hemisphere images required corrections of 0.775 mag- 
nitudes for both filters. This dependence on hemisphere 
may be related to errors in the rotational phase angle cor- 
rections. 

The next correction to the Voyager brightness values 
was for wavelength. The effective wavelengths of the clear 
and orange filters are 0.47 and 0.55 /urn, respectively. This 
corresponds to the wavelengths of the telescopic observa- 
tions. The blue and green Voyager filters have effective 
wavelengths of 0.48 and 0.58 gm, respectively. A wave- 
length correction was applied to the blue filter data to 
match it with the clear filter data by comparing the bright- 
ness values from these images at overlapping solar phase 
angles. A similar correction was applied to the green filter 
data to match it with the orange filter data. The filter 
corrections for the Europa phase curves were -0.104 and 
—0.074 magnitudes for the Voyager I blue filter images 
for the leading and trailing hemispheres, respectively. The 
corresponding Voyager II corrections were 0.076 and 0.092 
magnitudes. The filter corrections for the Europa phase 
curves for the orange filter Voyager I images were 0.015 
magnitudes for the trailing hemisphere, and no correction 
was applied to the leading hemisphere orange filter images. 
No corrections for wavelength were applied to the images 
used in the Ganymede phase curves. The corrections to 
the blue filter images for the Callisto phase curve measure- 
ments were 0.029 and 0.041 for the Voyager I and II images, 
respectively. Similarly for the Callisto orange filter images 
the corrections were 0.028 and —0.015 for Voyagers I and 
II, respectively. 

Once the camera corrections were applied then a similar 
comparison was made between the two spacecraft data 
sets for each filter. The Voyager I images were corrected 
to coincide with the Voyaer II images where necessary. 
These corrections consisted of —0.008 and —0.115 magni- 
tude corrections to the 0.55-^tm (orange and green) Europa 
phase curves for the leading and trailing hemispheres, re- 
spectively. The corresponding corrections for the 0.47 -/xm 
(clear and blue) Europa phase curves were 0.082 and 0.193, 
respectively. The correction to the Voyager I Ganymede 
images was 0.06 for both hemispheres and wavelengths. 
The Callisto phase curves required a correction of only 
0.125 to the Voyager II 0.47-/xm images to bring them into 
correspondence with the Voyager I images. 

The final correction to the Voyager data was made by 
comparing it to the corresponding telescopic observations. 
The clear and blue filter data were corrected to match the 
0.47-/xm telescopic data (in the case of Ganymede the 
telescopic observations were at 0.44 /xm) by comparing 
the brightness values at overlapping solar phase angles. 
Similarly, the orange and green filter data were corrected 
to match the 0.55-^m telescopic data. These corrections 
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for Europa's phase curves at 0.47 jum were -0.259 and 
-0.317 magnitudes for the leading and trailing hemi- 
spheres, respectively. The corresponding corrections to the 
Europa 0.55-/xm phase curves were —0.124 and —0.141 
magnitudes, respectively. For Ganymede the corrections 
were -0.12 and —0.14 magnitudes for both hemispheres 
at 0.47 and 0.55 ^ m, respectively. These corrections for 
Callisto’s phase curves at 0.47 ft m were —0.17 and —0.253 
magnitudes for the leading and trailing hemispheres, re- 
spectively, and were —0.07 and —0.157 magnitudes for 
Callisto’s phase curves at 0.55 fi m for the leading and 
trailing hemispheres, respectively. 

ANALYTICAL METHODS 

In the first set of analyses concerning the affect of the 
choice of single particle scattering function, Hapke’s model 
(Hapke 1981, 1984, 1986) was fit to the data using a modi- 
fied least squares grid search. Each data set was modeled 
twice, once using Hapke’s model with a two parameter 
double Henyey-Greenstein function (2P-HG function) 
and a second time using a three parameter double Henyey- 
Greenstein function (3P-HG function). Isotropic multiple 
scattering was assumed. The form of the 2P-HG function 
used is given by 

P 2 (a) = ((1 - c 2 )( 1 - bl)![ 1 + 2b 2 cos (a) + b\] y2 ) 

+ (c 2 ( l — b 2 )/[\ ~ 2fr 2 cos (a) + b 2 ] 312 ), 

where a is the phase angle, and b 2 and c 2 are the single 
particle scattering parameters. The form of the 3P-HG 
function used is given by 

P,(a) = ((1 - c 3 )(l - hi )/ (1 T 2b, cos (a) + b 2 ,f 2 ) 

4- (c,( 1 — d 2 )/[\ — 2d, cos (a) T dl] V2 ), 

where a: is the phase angle, and b 3 , c 3 , and d, are the 
single particle scattering parameters. Note that the 3P-HG 
function defaults to the 2P-HG function for cases, where 
d, = ~b,. The results of fitting Hapke’s isotropic multiple 
scattering model to our data sets are listed in Table II along 
with the corresponding RMS values. The RMS values are 
defined such that 


RMS 


Sw,- m„y 



where N is the number of observational data points, M is 
the calculated magnitude from Hapke’s model, m is the 
observed magnitude, and the sum is over the number of 
points. Figures 2 through 7 plot our theoretical results 
versus the data sets. 

In the second set of analyses Hapke’s model was fit to 


the phase curves of Europa and Ganymede using the 
method of Verbiscer (1991), which incorporates aniso- 
tropic multiple scattering. This method assumes that the 
anisotropy in the single and multiple scattered components 
is described by the same scattering function. The 3P-HG 
function was used. The results are listed in Table III and the 
corresponding phase curves are plotted in Figs. 2 through 7. 
Only the 0.47-/zm Ganymede data were fit since it has 
adequate large (>90°) phase angle coverage. The Callisto 
phase curves were not modeled since its surface is darker 
than Europa’s or Canymede’s; thus the role of multiple 
scattering is weaker. Anisotropic multiple scattering will 
be most important for bright surfaces, therefore any differ- 
ences in isotropic versus anisotropic multiple scattering 
will be most evident in the modeling of the brighter surfaces 
of Europa and Ganymede. 

ERROR ANALYSIS 

Error analysis based on the data set and on the mathe- 
matical inter-relationship of the Hapke model parameters 
is complicated and non-trivial. The value, and changes in 
the value, of each parameter affects the shape of the phase 
curve in different phase angle ranges. Changes in the shape 
of the phase curve due to changes in the value of one 
parameter oftentimes may be countered by changes in 
another parameter. Therefore, while an error analysis be- 
gins with studying the affects of the individual parameters 
on the shape of the phase curve, the more convoluted 
variations in phase curve shape due to changes in multiple 
parameters must eventually be considered and evaluated. 
Below is a description of the evaluation of the error bars 
on the Hapke parameter values for each satellite. These 
error bars are summarized in Table IV. 

A first order analysis was accomplished by comparing 
the phase curve of a standard set of Hapke parameters 
to phase curves where each parameter has been varied 
individually while the others have been kept constant in 
order to observe the variations in the phase curve due to 
the change in that one parameter. It can be shown that 
the value of the single scattering albedo, w, affects the 
entire range of phase angles, whereas the opposition pa- 
rameters bo (amplitude) and h (width) affect the shape of 
the phase curve from 0° to 30° solar phase angle. The 
surface roughness parameter, 0- bar, affects the shape of 
the phase curve most dramatically at solar phase angles 
greater than 90°; however, it can influence the shape of 
the phase curve at phase angles as small as 20°. The affects 
of the single particle scattering function also span the entire 
range of phase angles where the values of the single scatter- 
ing function parameters affect different ranges of phase 
angle depending on the scattering function chosen. Initial 
error bars are set by comparing the variation between two 
phase curves, where the value of one parameter has been 
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TABLE II 

Hapke Parameters for the Icy Galilean Satellites and Rhea 


Leading Trailing Leading Trailing 



0.47 jiun 

0.55 fim 

0.47 fim 

0.55 ptm 


0.47 fim 

0.55 ijl m 

0.47 fim 

0.55 fim 

w 

0.922 

0.964 

0.897 

Europa 

0.940 

w 

0.934 

0.964 

0.897 

0.930 

bo 

0.50 

0.50 

0.45 

0.50 

bo 

0.49 

0.43 

0.51 

0.521 

h 

0.0016 

0.0016 

0.0016 

0.0016 

h 

0.0015 

0.0016 

0.0016 

0.0016 

b2 

0.431 

0.429 

0.43 

0.443 

b3 

0.770 

0.726 

0. 

0.083 

c2 

0.921 

0.887 

0.713 

0.609 

c3 

0.780 

0.945 

0.691 

0.784 

0-bar 

10 

10 

10 

10 

d3 

0-bar 

-0.459 

8 

0.416 

10 

-0.417 

11 

-0.386 

11 

rms 

0.029 

0.029 

0.039 

0.033 

rms 

0.029 

0.029 

0.035 

0.025 

w 

0.830 

0.945 

0.870 

Ganymede 

0.81 

w 

0.830 

0.930 

0.870 

0.810 

bo 

0.62 

0.86 

LOO 

0.23 

bo 

0.63 

0.91 

LOO 

0.230 

h 

0.003 

0.004 

0.074 

0.074 

h 

0.003 

0.003 

0.074 

0.074 

b2 

0.282 

0.380 

0.039 

0.307 

b3 

0.090 

0.200 

0 

0.350 

c2 

0.960 

0.427 

0.989 

0.962 

c3 

0.960 

0.450 

0.820 

0.970 

tf-bar 

28 

29 

35 

35 

d3 

0-bar 

-0.280 

28 

- 0.380 
28 

-0.046 

35 

-0.305 

35 

rms 

0.067 

0.036 

0.097 

0.047 

rms 

0.067 

0.036 

0.097 

0.047 

w 

0.740 

0.540 

0.470 

Callisto 

0.550 

vv 

0.510 

0.605 

0.530 

0.650 

bo 

1.00 

1 .00 

0.27 

0.73 

bo 

0.82 

1 .00 

0.50 

0.87 

h 

0.031 

0.031 

0.0277 

0.0277 

h 

0.031 

0.031 

0.0277 

0.0277 

b2 

0.729 

0.132 

0.432 

0.206 

b3 

0 

0.005 

0.787 

0.988 

c2 

0.024 

0.949 

0.542 

0.958 

c3 

0.033 

0.015 

0.489 

0.618 

0-bar 

42 

42 

42 

42 

d3 

0-bar 

-0.694 

42 

0.687 

42 

-0.400 

42 

-0.238 

42 

rms 

0.093 

0.062 

0.213 

0.066 

rms 

0.092 

0.061 

0.195 

0.065 

Note. 

Error bars: vv = 

±0,01. bo = ±0.01, h 

= ±0.001 , b2 and b3 = ±0.005, 

c2 and c3 - 

±0.005. d3 - 

- ±0.005, fl-bar = 

±2. 



varied, with the scatter and phase angle coverage in the 
data sets. 

The initial error bars for the parameters w\ bo. and h 
are ±0.02, ±0.1, ±0.002, respectively, for all three of the 
icy Galilean satellites. The initial error bars for fl-bar are 
determined by the presence of large phase angle data. If 
data exist beyond phase angles of 100° the initial error 
bars are ±5°. If the large phase angle data exists only out 
to 70° the initial error bars are ±10°, and similarly, if the 
phase angle coverage exists only out to 50° the initial error 
bars are ±15°. The Europa the initial error bars in 0-bar 
are ±5° for the leading hemisphere in both wavelengths 
and for the trailing hemisphere at 0.47 /xm. Using these 
criteria the initial error bar for Europa’s trailing hemi- 
sphere at 0.55 fim is ±10°. However, the surface roughness 
parameter should be independent of wavelength; therefore 
the error bar on this parameter for Europa’s trailing hemi- 
sphere at 0.55 fim is also constrained by the modeling 
results to the 0.47-/xm data for this hemisphere. Therefore 
the error bar is less than ±10° and closer to ±5°. The initial 
error bars for the 0.47-jum data for both hemispheres of 
Ganymede are ±5°, since there is a data point at 120° 


phase angle and the scatter in the large phase angle mea- 
surements of the other satellites correlate to 5° of variaiton 
in tf-bar. For the 0.55-/xm Ganymede data sets the initial 
error bars based on the range of large phase angle data 
are ±15°. However, since tf-bar is independent of wave- 
length, the modeling results of the 0.47-/xm data constrain 
the modeling solutions to the 0.55-/xm data (and vice 
versa). Therefore the error bars for the Ganymede 0.55- 
/xm data are closer to ±10° than ±15°. The initial error 
bars for 0-bar for Callisto are ±5° for the trailing hemi- 
sphere data at both wavelengths and ±15° for the leading 
hemisphere in both wavelengths, based on the phase angle 
coverage in the data sets. 

In the analyses presented in this paper two different 
single particle scattering functions were used; the two pa- 
rameter double Henyey-Greenstein function (2P-HG), 
with parameters b 2 and r>, and the three parameter double 
Henyey-Greenstein function (3P-HG), with parameters 
c 3 , and d 3 . The 2P-HG function assumes that the half- 
width of the forward and back scattering lobes are equiva- 
lent and that the relative amplitude varies. The parameter 
b 2 is the width parameter and the parameter c 2 is the 
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Europa Leading Hemisphere at 0.47 microns Europa Leading Hemisphere at 0.47 microns 






phase ongle (degrees) phase angle (degrees) 

FIG. 2. The solar phase curves of Europa's leading (a and b) and trailing (c and d) hemispheres at 0.47 fx m, where the diamonds represent 
the disk-integrated observations, the solid line represents the isotropic multiple scattering Hapke model solution incorporation a two parameter 
HG function, the dashed line represents the isotropic multiple scattering Hapke model solution incorporating a three parameter HG function, and 
the dash-dot line represents the anisotropic multiple scattering Hapke model solution incorporation a three parameter HG function. Those places 
where only a solid line is seen indicates that the solutions lie on top of each other and, therefore, cannot be distinguished. 


amplitude parameter. Variations in each of these parame- 
ters affects the shape of the phase curve at all phase angles. 
Any changes at large phase angles, where there are fewer 
data to constrain the phase curves studied here, also result 
in changes at low phase angles which is well constrained 
by the data sets in this study. While changes in the shape 
of the phase curve due to variations in the values of b 2 
and c 2 can often by negated by changes in the values of 
0-bar and the opposition parameters, for a given set of b 2 
and c 2 there is a correlating set of unique values for 0-bar 
and the opposition parameters (and vice versa). Using the 
method described above the initial error bars for b 2 and 
c 2 for each of the icy Galilean satellites are ±0.02. 

The second single particle scattering function used in this 
study is the three parameter double Henyey-Greenstein 
function (3P-HG), where b 3 is the half-width of the forward 
scattering component, d 3 is the half-width of the backscat- 


tered component, and c 3 is the relative amplitude of these 
two scattering components. The values of b 3 are con- 
strained by the high phase angle (a > 90°) data and the 
values of d 3 are constrained by the low phase angle (a < 
90°) data. The value of c 3 affects the shape of the phase 
curve over the entire range of phase angles. Changes at 
large phase angles due to variations in the values of b 3 or 
c 3 can oftentimes be negated by changes in the value of 
0-bar depending on the magnitude of the variations in the 
scattering function parameters. Changes in the shape of 
the phase curve at small phase angles due to changes in 
the values of d 3 or c 3 can sometimes be compensated for 
by changes in the opposition parameters bo and h . Using 
the same methods as above the initial error bars for d 3 , 
and c 3 are ±0.01 and ±0.02, respectively, for the all three 
of the icy Galilean satellites. The initial error bars for b 3 
depend on the presence and extent of the large phase angle 
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Europa Leading Hemisphere at 0.55 microns 




Europa Trailing Hemisphere at 0.55 microns 



phase angle (degrees) 
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FIG. 3. The solar phase curves of Europa’s leading (a and b) and trailing (c and d) hemispheres at i).55 /xm. where the diamonds represent 
the disk-integrated observations, the solid line represents the isotropic multiple scattering Hapke model solution incorporating a two parameter 
HG function, the dashed line represents the isotropic multiple scattering Hapke model solution incorporating a three parameter HG function, and 
the dash-dot line represents the anisotropic multiple scattering Hapke model solution incorporating a three parameter HG function. Those places 
where only a solid line is seen indicates that the solutions lie on top of each other and, therefore, cannot be distinguished. 


data. While the values for tf-bar also depend on large phase 
angle data, 0-bar has the added constraint that it is indepen- 
dent of wavelength. The parameter by , however, does vary 
with wavelength. Therefore the initial error bars for by 
are ±0.05 for the leading hemisphere of Europa at both 
wavelengths, the trailing hemisphere of Europa at 0.47 
ptm, both hemispheres of Ganymede at 0.47 jwm, and the 
trailing hemisphere of Callisto at both wavelengths. The 
initial error bars for by for Europa’s trailing hemisphere 
at 0.55 fim is ±0.1 and is ±0.5 for both hemispheres of 
Ganymede at 0.55 /xm and Callisto’s leading hemisphere 
at both wavelengths. 

Further constraints on the error bars of the Hapke pa- 
rameters can be placed by understanding the physical na- 
ture of the parameters. The parameter ft-bar is a measure 
of surface roughness (the tilt of the surface averaged over 
resolutions ranging from several particle sizes to the resolu- 
tion of the detector), which is a physical characteristic 


independent of wavelength. By analyzing the phase curves 
of the same surface at different wavelengths this parameter 
can be better constrained and its value can be decoupled 
from the values of the other parameters in Hapke's model. 
The value of 0-bar is mathematically coupled to the values 
of w and the single particle scattering function. Any 
method of analysis which decouples their affect on the 
shape of the phase curve also affects the error analysis. 
The comparison of the disk-integrated solutions to the 
parameters based on the disk-integrated phase curves to 
disk-resolved measurements of the same surface is an addi- 
tional method for decoupling the value of tf-bar, h\ and 
the single particle scattering function parameters. By using 
the disk-integrated parameter solutions and applying their 
values to disk-resolved data and using the wavelength inde- 
pendent nature of 0-bar, the error bars on w can be con- 
strained to ±0.01 based on the scatter in the disk-resolved 
data for all the icy Galilean satellites. 
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Ganymede Leading Hemisphere at 0.44 microns 



Ganymede Leading Hemisphere at 0.44 microns 



Ganymede Trailing Hemisphere at 0.44 microns 



Ganymede Trailing Hemisphere at 0.44 microns 



FIG. 4. The solar phase curves of Ganymede’s leading (a and b) and trailing (c and d) hemispheres at 0.47 />tm. where the diamonds represent 
the disk-integrated observations, the solid line represents the isotropic multiple scattering Hapke model solution incorporating a two parameter 
HG function, the dashed line represents the isotropic multiple scattering Hapke model solution incorporating a three parameter HG function, and 
the dash-dot line represents the anisotropic multiple scattering Hapke model solution incorporating a three parameter HG function. Those places 
where only a solid line is seen indicates that the solutions lie on top of each other and, therefore, cannot be distinguished. 


The opposition parameter h is dependent on the grain 
size distribution and porosity within the optically active 
portion of the regolith. These physical characteristics are 
also independent of the wavelength of observation. Analy- 
sis of the phase curves of the same surface at different 
wavelengths in addition to good quality disk-integrated 
measurements at very low phase angles ( a < 1°) constrains 
the values of this parameter to wihin ±0.0005 for all the icy 
satellite data except for Ganymede’s trailing hemisphere 
which lacks the very low phase angle data. The value of 
the opposition width parameter, /z, is strongly coupled to 
the value of bo, the opposition amplitude. Using the wave- 
length independent nature of h and the quality of the very 
low phase angle data, the value of bo can be constrained 
to ±0.02 for all the data sets excluding Ganymede’s trail- 
ing hemisphere. 

Additional constraints on the error bars for the Hapke 


parameters may also be based on the method of analysis. 
For the data sets fits using a modified least squares grid 
search, differences in RMS values of 1% are considered 
significant and error bars can be established on this basis. 
The modified least squares grid search is an iterative pro- 
cess where the initial grid size is coarse and the solutions 
with the 10 least RMS values are all considered viable 
solutions to the model. A second grid search is run where 
the limits of the grid for a specific parameter are based on 
the range of values found from the first 10 solutions plus 
and minus the initial grid increment and the grid size is 
reduced. Once again the solutions with the 10 least RMS 
values are considered viable solutions to the model and 
the process is repeated until either the desired grid size is 
reached or the differences in the least 10 RMS values is 
less than 1%. For example, the parameter w will initially 
be run from values of 0.2 to 1.0 in increments of 0.2. The 
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Ganymede Leading Hemisphere at 0.55 microns 
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FIG. 5. The solar phase curves of Ganymede’s leading (a and b) and trailing (c and d) hemispheres at 0.55 jum, where the diamonds represent 
the disk-integrated observations, the solid line represents the isotropic multiple scattering Hapke model solution incorporating a two parameter 
HG function, and the dashed line represents the isotropic multiple scattering Hapke model solution incorporating a three parameter HG function. 
Those places where only a solid line is seen indicates that the solution lie on top of each other, and therefore, cannot be distinguished. 


range of values for w from the first set of 10 lesat RMS 
solutions may range from 0.6 to 0.8. Then the next grid 
search run will have values for w ranging from 0.4 to t.O 
in increments of 0.1. All the parameters are treated in 
the same manner. The error bars listed in Table II 
correspond to the smallest grid increment run on that pa- 
rameter. However, by examining the parameter values for 
the 10 least RMS solutions along each step in the grid 
search, coupled with the error analysis results discussed 
above, the following final error bars can be set for the 
Hapke parameter values (in the isotropic multiple scatter- 
ing case) for the icy Galilean satellites: w = ±0.01 for 
all the satellites, bo = ±0.02 for all the satellites except 
Ganymede's trailing hemisphere where the value is ±0.1, 
h = ±0.0005 for all the satellites except Ganymede’s trail- 
ing hemisphere where the value is ±0.01, 0-bar = ±5° 
except for Ganymede’s 0.55-^m phase curves where the 
error bars are ±10° and Callisto’s leading hemisphere 


phase curves where the tf-bar error bars are ±15°. The 
final error bars for the scattering function parameters are 
b 2 * C2 = C3 = ±0.02 and d3 = ±0.01. The error bars for 
by are ±0.05 for Europa’s OAl-jxm phase curves, Europa’s 
0.55-ftm leading hemisphere phase curve, Ganymede’s 
0.47 -/lhh phase curves, and Callisto’s trailing hemisphere 
phase curves. For Europa’s 0.55 -/xm trailing hemisphere 
phase curve the by error bar is ±0.1 and is ±0.5 for Ga- 
nymede’s 0.55-/xm phase curves and Callisto’s leading 
hemisphere phase curves. 

The phase curve analyses incorporating anisotropic mul- 
tiple scattering were conducted using software designed 
by Paul Helfenstein following the method outlined by Ver- 
biscer (1991). The error analysis in the software utilized a 
gradient-expansion search of the \ 2 parameter space which 
locates the minimum between measured and calculated 
reflectance values. Details of this method are given in Be- 
vington (1969) and Helfenstein (1986). 
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Callisto Leading Hemisphere at 0.47 microns 




0 50 100 150 

phase angle (degrees) 


Caifisto Leading Hemisphere at 0.47 microns 



Cailisto Trailing Hemisphere at 0.47 microns 



phase angle (degrees) 


FIG. 6. The solar phase curves of Callisto’s leading (a and b) and trailing (c and d) hemispheres at 0.47 /zm, where the diamonds represent 
the disk-integrated observations, the solid line represents the isotropic multiple scattering Hapke model solution incorporating a two parameter 
HG function, the dashed line represents the isotropic multiple scattering Hapke model solution incorporating a three parameter HG function. Those 
places where only a solid line is seen indicates that the solutions lie on top of each other and, therefore, cannot be distinguished. 


RESULTS 

Isotropic Multiple Scattering Results 

General. Hapke's theory models the phase curves of 
these satellites very well, as demonstrated in Figs. 2 through 
7. There is no substantial difference in the quality of fit 
between the model incorporating the 3P-HG function over 
the model incorporating the 2P-HG function. The largest 
difference in RMS values is less than 2%, where as the 
majority of fits are within 0.1 %. There were only four phase 
curves for which no large phase angle ( a > 60°) data were 
available. In the case of the 0.55-pm data for both hemi- 
spheres of Ganymede the two models predict the same 
behavior at large phase angles. However, for the leading 
hemisphere phase curves of Callisto the two models predic- 
tions at large phase angles diverge significantly. Therefore 
the results of the models for these phase curves must be 
interpreted with caution. 


Figures 8 through 10 are comparisons of both model 
predictions versus disk-resolved scans across the satellite 
surfaces. Each scan is a vertical scan from north to south 
across an orthographic projection of the image where the 
values of incidence, emission, latitude, and longitude 
angles and pixel brightness are recorded based on the im- 
age label information. Figure 8 is a scan across Europa 
where both models fit the lower albedo data (which corre- 
sponds to mottled terrains), but do not model the bright 
plains units at low and high pixel locations. This type of 
fit is to be expected when albedo varies across the surface. 
Figure 9 is a scan across Ganymede where the two models 
made very different predictions of the surface brightness 
as a function of location, but each fits the data equally 
well. The 2P-HG model fits the higher albedo regions, 
whereas the 3P-HG model fits the lower albedo regions. 
Once again, these types of fits are to be expected when 
albedo varies across the surface. Figure 10 is a scan across 
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Callisto Leading Hemisphere at 0.55 microns 
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Callisto Trailing Hemisphere at 0.55 microns 
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Callisto Leaaing Hemisphere at 0.55 microns 



FIG. 7. The solar phase curves of Callistcfs leading (a and b) and trailing (c and d) hemispheres at 0.55 pirn, where the diamonds represent 
the disk-integrated observations, the solid line represents the isotropic multiple scattering Hapke model solution incorporating a two parameter 
HG function, the dashed line represents the isotropic multiple scattering Hapke model solution incorporating a three parameter HG function. Those 
places where only a solid line is seen indicates that the solutions lie on top of each other and, therefore, cannot be distinguished. 


the surface of Callisto where both models fit the data 
equally well. 

Roughness. The macroscopic roughness parameter, 6- 
bar, shows no hemispherical differences on Europa or Cal- 
listo in both models. This agrees with the results Domingue 
etal (1991) obtained for Europa, but contradicts the results 
Buratti (1991) obtained for Callisto. Buratti (1991) found 
that the leading hemisphere of Callisto is macroscopically 
rougher than its trailing hemisphere. The discrepancy be- 
tween our results and those of Buratti (1991) is due to 
the fact that neither study had large phase angle data for 
Callisto’s leading hemisphere. In Hapke’s model it is essen- 
tial to have large phase angle data in order to constrain 
the macroscopic roughness parameter. Since there are cur- 
rently no large phase angle data for Callisto’s leading hemi- 
sphere, it cannot be stated with certainty if there is a hemi- 
spherical difference in surface roughness on Callisto. The 
shape of the single scattering function is also greatly influ- 


enced by large phase angle observations. Therefore the 
different choices in single scattering function coupled with 
the lack of large phase angle data also contribute to the 
different roughness results between the two studies. No 
differences in the value of roughness were seen for Callisto 
between the 2P-HG model and the 3P-HG model. The 
values of 0-bar for Europa between the two models are 
within 1° of each other, well within the error bars. Both 
models show that the trailing hemisphere of Ganymede is 
rougher than the leading hemisphere by about 6°. Buratti 
(1991) found no hemispherical differences in roughness in 
her analysis of Ganymede, which we attribute to the choice 
of single particle scattering function. Hartman and Dom- 
ingue (1997) found that the two models returned very 
different values of 0-bar when applied to the lunar phase 
curve, up to 12° difference. However, no differences within 
each satellite were found outside the error bars. An in- 
creasing trend in the 0-bar values with distance from Jupi- 
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TABLE III 

Anisotropic Multiple Scattering Hapke Model Results 




Leading 



Trailing 


0.47 fi 

Error bar 

0.55 fxm 

Error bar 

0.47 jiim 

Error bar 

0.55 /txm 

Error bar 

w 

0.905 

±0.001 

0.949 

±0.002 

0.897 

±0.024 

0.911 

±0.002 

bo 

0.45 

±0.19 

0.5 

±0.17 

0.74 

±0.19 

0.301 

±0.041 

h 

0.0019 

±0.006 

0.0013 

±0.0029 

0.0009 

±0.0015 

0.02 

±0.005 

b3 

0.5 

nc 

0.0 

±0.1 

0.0 

nc 

0.0003 

±0.01 

c3 

0.972 

±0.017 

0.967 

±0.04 

0.681 

±0.102 

0.988 

±0.007 

d3 

-0.43 

±0.01 

-0.4185 

±0.01 

-0.418 

±0.107 

-0.333 

±0.001 

fl-bar 

2. 

±5 

0.1 

±5 

13.7 

±2.6 

0 

±5 

rms 

0.027 


0.03 


0.035 


0.023 






Ganymede 





w 

0.8227 

±0.02 

n/a 

n/a 

0.873 

±0.016 

n/a 

n/a 

bo 

0.5537 

±0.1762 

n/a 

n/a 

1.2154 

±0.2722 

n/a 

n/a 

h 

0.0035 

±0.0045 

n/a 

n/a 

0.0758 

±0.016 

n/a 

n/a 

b3 

0.01 

nc 

n/a 

n/a 

0.0395 

nc 

n/a 

n/a 

c3 

-0.92 

±0.1 

n/a 

n/a 

0.5592 

±0.12 

n/a 

n/a 

d3 

-0.29 

±0.09 

n/a 

n/a 

-0.0422 

±0.05 

n/a 

n/a 

fl-bar 

27.4 

±2.6 

n/a 

n/a 

35.06 

±1.6 

n/a 

n/a 

rms 

0.067 




0.096 





Note . nc, not calculated, na, not available. 


ter is seen, which correlates with the amount of large im- 
pact cratering observed on these satellites and therefore 
the age of their surfaces. 

Opposition effect. The opposition width parameter, /?, 
showed no differences outside the error bars between the 
2P-HG model and the 3P-HG model, consistent with the 
findings of Hartman and Domigue (1997), No hemispheri- 
cal differences were seen on Europa; however, both Ga- 
nymede and Callisto displayed variations with hemisphere. 
The value of the opposition width parameter is dependent 
on the porosity and grain size distribution of the particles 
within the optically active portion of the regolith. Gener- 
ally the larger the value of this parameter the more compact 


or less porous the regolith. The correlation between h and 
grain size is less straightforward. Therefore the values of 
h seen for both hemispheres of Europa indicate a uniform, 
highly porous, surface. The results shown in Table II indi- 
cate that the leading hemisphere of Ganymede is much 
more porous than the trailing hemisphere. The surprising 
result is that the value of h for the trailing hemisphere of 
Ganymede is even larger than that obtained by Hartman 
and Domingue (1997) for the Moon! Buratti (1991) found 
no hemispherical differences in the compaction state of 
Ganymede’s regolith. The results for Callisto show that its 
leading hemisphere is slightly more compact than its trail- 
ing hemisphere. Buratti (1991) found the opposite; how- 


Satellite 


Europa 


Ganymede 


TABLE IV 

Isotropic Model Error Bars 


Hemisphere 

w 

bo 

h 

b 2 

c 2 

b 3 

d 3 

C 3 

0-bar 

0.47-/mm lead 

0.01 

0.02 

0.0005 

0.02 

0.02 

0.05 

0.01 

0.02 

5 

0.47-zxm trail 

0.01 

0.02 

0.0005 

0.02 

0.02 

0.05 

0.01 

0.02 

5 

0.55-/xm lead 

0.01 

0.02 

0.0005 

0.02 

0.02 

0.05 

0.01 

0.02 

5 

0.55ju.m trail 

0.01 

0.02 

0.0005 

0.02 

0.02 

0.1 

0.01 

0.02 

5 

0.47-^tm lead 

0.01 

0.02 

0.0005 

0.02 

0.02 

0.05 

0.01 

0.02 

5 

0.47-^m trail 

0.01 

0.1 

0.01 

0.02 

0.02 

0.05 

0.01 

0.02 

5 

0.55-ftm lead 

0.01 

0.02 

0.0005 

0.02 

0.02 

0.5 

0,01 

0.02 

10 

0.55/xm trail 

0.01 

0.1 

0.01 

0.02 

0.02 

0.5 

0.01 

0.02 

10 

0.47-ptm lead 

0.01 

0.02 

0.0005 

0.02 

0.02 

0.5 

0.01 

0.02 

15 

0.47-^tm trail 

0.01 

0.02 

0.0005 

0.02 

0.02 

0.05 

0.01 

0.02 

5 

0.55-/xm lead 

0.01 

0.02 

0.0005 

0.02 

0.02 

0.5 

0.01 

0.02 

15 

Q.55fim trail 

0.01 

0.02 

0.0005 

0.02 

0.02 

0.05 

0.01 

0.02 

5 


Callisto 
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FIG. 8. Vertical scan from north to south across an orthographic projection of the Voyager image 16342.22 of Europa. This image is an orange 
filter image of the trailing hemisphere. The diamonds represent the Voyager data, the solid line represents the isotropic multiple scattering 2P-HG 
Hapke model solution, and the dashed line represents the isotropic multiple scatteing 3P-HG Hapke model solution. 


ever, our phase curve contains opposition data which was 
not available at the time she performed her study. The 
differences in the results obtained by Buratti (1991) and 
this study for Ganymede are not understood and may be 


attributed to the lack of the very low phase angle data 
needed to well define the opposition effect. 

There are variations in the value of the opposition ampli- 
tude parameter, bo, between the 2P-HG function and the 



FIG. 9. Vertical scan from north to south across an orthographic projection of the Voyager image 20635.35 of Granymede. This image is an 
orange filter image of the trailing hemisphere. The diamonds represent the Voyager data, the solid line represents the isotropic multiple scattering 
2P-HG Hapke model solution, and the dashed line represents the isotropic multiple scattering 3P-HG Hapke model solution. 
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FIG. 10. Vertical scan from north to south across an orthographic projection of the Voyager image 20616.1 1 of Callisto. This image is an orange 
filter image of the trailing hemisphere. The diamonds represent the Voyager data, the solid line represents the isotropic multiple scattering 2P-HG 
Hapke model solution, and the dashed line represents the isotropic multiple scattering 3P-HG Hapke model solution. 


3P-HG function, which was not found by Hartman and 
Domingue (1997) in their similar analysis of the Moon. 
In some cases the difference is substantial, in others the 
differences are within the error bars of the parameters. 
The opposition parameter bo is correlated to the opacity 
of the regolith particles. It is a measure of the percentage 
of incident light which is scattered at the surface of the 
particle compared to the total amount of light scattered 
(both from the surface and internally). A bo value of one 
indicates an opaque particle where all the light is scattered 
at the particle surface. In the case of Europa both models 
return values of bo near one-half for both hemispheres, 
indicating a fairly transparent material, which is consistent 
with a predominately water ice surface. The results for 
Ganymede show variations between hemispheres and be- 
tween wavelengths. Both models indicate that on the lead- 
ing hemisphere the material is more opaque at 0.55 gm 
than at 0.44 gm, while the opposite is true on the trailing 
hemisphere. The trailing hemisphere values of bo do not 
vary between the two models; however, there is a signifi- 
cant difference in the bo values returned by the models 
for the leading hemisphere at 0.55 gm. The 3P-HG model 
indicates a more opaque material than the 2P-HG model. 
The results for Callisto also demonstrate variations be- 
tween hemispheres and wavelengths. Both models show 
that the material on the leading hemisphere is more opaque 
than the material on the trailing hemisphere; however, the 
2P-HG model indicates larger differences in opacity at 0.47 
gm than the 3P-HG model. The 2P-HG model indicates 


significantly more opaque material at 0.47 /xm on the lead- 
ing hemisphere than the 3P-HG model, while it indicates 
significantly more transparent material at both wave- 
lengths on the trailing hemisphere than the 3P-HG model. 

Single scattering parameters. Hapke 's theory also re- 
turns a measure of the single scattering albedo, w. The 
values of vv vary significantly between the two model re- 
sults. The 2P-HG model finds a lower w value for Europa’s 
leading hemisphere at 0.47 /xm, but returns a similar value 
as the 3P-HG model for Europa’s leading hemisphere at 
0.55 /xm. The 2P-HG model also finds similar values as 
the 3P-HG model for Europa’s trailing hemisphere at 0.47 
/xm, but returns a brighter value for this hemisphere at 
0.55 /xm. The only differences in vv found by these two 
models for the hemispheres of Ganymede is for the leading 
hemisphere at 0.55 /xm, in which the 2P-HG model finds 
a brighter value than the 3P-HG model. The differences for 
Callisto are very large. The 3P-HG model finds significantly 
brighter vv values at both wavelengths for the trailing hemi- 
sphere and at 0.55 /xm for the leading hemisphere. How- 
ever the 3P-HG model finds a darker value for the leading 
hemisphere at 0.47 /xm than the 2P-HG model. 

Figures 11 through 13 plot and compare the single scat- 
tering functions from the models for each satellite. Both 
isotropic multiple scattering models show that Europa’s 
leading hemisphere is backscattering at both wavelengths; 
however, the isotropic 3P-HG model also predicts a sig- 
nificant forward scattered component. In the case of Euro- 
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FIG. 11 . The single particle scattering functions for Europa’s (a) leading and (b) trailing hemispheres plotted versus phase angle. The long and 
short dashed lines represent the isotropic multiple scattering model with the two parameter HG function solutions at 0.47 and 0.55 gm. respectively. 
The dash-dot and dash-double-dot lines represent the isotropic multiple scattering model with the three parameter HG function solutions at 0.47 
and 0.55 gm. respectively. The thin and thick solid lines represent the anisotropic multiple scattering model with the three parameter HG function 
solutions at 0.47 and 0.55 gim, respectively. Those areas where only one solution is apparent indicate areas where the scattering functions lie on 
top of each other and, therefore, cannot be disginguished. 


pa’s trailing hemisphere the isotropic 2P-HG function 
shows a significant forward scattering lobe. The isotropic 
3P-HG function, on the other hand, shows no evidence of 
forward scattering on Europa’s trailing hemisphere in ei- 
ther wavelength. Both isotropic scattering models show 
that both hemispheres of Europa have a significant back- 


scattering component, with the amplitude of the backscat- 
tering lobe stronge on the leading hemisphere. 

Figure 12 shows the single particle scattering behavior 
for Ganymede. Both isotropic 2P-HG and 3P-HG func- 
tions show some interesting hemispherical and wavelength 
dependent differences. At 0.47 fim the single scattering 
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FIG. 12. The single particle scattering functions for Ganymede’s (a) leading and (b) trailing hemispheres plotted versus phase angle. The long 
and short dashed lines represent the isotropic multiple scattering model with the two parameter HG function solutions at 0.47 and 0.55 ftm, 
respectively. The dash-dot and dash-double-dot lines represent the isotropic multiple scattering model with the three parameter HG function 
solutions at 0.47 and 0.55 /u m, respectively. The solid lines represent the anisotropic multiple scattering model with the three parameter HG function 
solutions at 0.47 /xm. Those areas where only one solution is apparent indicate areas where the scattering functions lie on top of each other and, 
therefore, cannot be distinguished. 


function predicted by both models is isotropic on the trail- 
ing side, while the leading hemisphere it is backscattering. 
At 0.55 fi m on the trailing hemisphere both isotropic multi- 
ple scattering models show a backscattering scattering 
function; however, for the leading hemisphere the isotropic 
2P-HG model predicts a slightly stronger forward scatter- 


ing component. The isotropic 3P-HG function at 0.55 /xm 
predicts that both hemispheres are predominately back- 
scattering with little difference between the two. However, 
there was no large phase angle data for the leading hemi- 
sphere of Ganymede at 0.55 /xm, so the results for the 
leading hemisphere are preliminary. These findings dis- 
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FIG. 13. The single particle scattering functions for Callisto’s (a) leading and (b) trailing hemispheres plotted versus phase angle. The solid 
and dashed lines represent the isotropic multiple scattering model with the two parameter HG tunetion solutions at 0.47 and 0.55 ^m. respectively. 
The dash-double-dot and dash-dot lines represent the isotropic multiple scattering model with the three parameter HG function solutions at 0.47 
and 0.55 jim, respectively. Those areas where only one solution is apparent indicates areas where the scattering functions lie on top ol each other 
and, therefore, cannot be distinguished. 


agree with Burattfs (1991) results that there are no hemi- 
spherical differnces in the single scattering behavior on 
Ganymede; however, this may be an artifact of the different 
single particle scattering functions chosen between the 
two studies. 

Figure 13 shows the single particle scattering behavior 
for Callisto. Recall that there were no large phase angle 


data for the leading hemisphere at either wavelength, so 
the leading hemisphere results need to be interpreted with 
great caution. The isotropic multiple scattering 3P-HG 
function predicts that the leading hemisphere’s particles 
are nearly isotropic, whereas the trailing hemisphere is 
more backscattering with some forward scattering in the 
0.47-/xm phase curve results. The isotropic 2P-HG function 
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predicts a highly forward scattering surface for both hemi- 
spheres at both wavelengths. This large difference is the 
result of no large phase angle measurements in the solar 
phase curve. Buratti (1991) found that the leading hemi- 
sphere of Callisto was more backscattering than the trailing 
hemisphere; however, her data set also did not include 
large phase angle observations for the leading hemisphere. 
Thus no conclusive statement can be made regarding hemi- 
spherical differences in the single scattering behavior on 
Callisto. 

When interpreting the single particle scattering results 
a few caveats must be stated. McGuire and Hapke (1995) 
and Hartman and Domingue (1997) both demonstrated 
that the single scattering functions used in the analyses 
presented here are viable functions and well describe the 
laboratory scattering data of McGuire and Hapke. There- 
fore the significance between the isotropic 2P-HG and 3P- 
HG results need to be explained. 

The effects of a backward scattering component to the 
scattering function are most strongly seen between 0 and 
40° phase angle. Likewise, the effects of a forward scatter- 
ing component to the scattering function are most strongly 
seen beyond 140° phase angle. Since none of the phase 
curves analyzed here has measurements beyond 140°, the 
forward scattering component of the scattering function is 
not well constrained (see Error Analysis). However, these 
two scattering functions do not treat scattering in the same 
manner. The half-width of the forward and backward scat- 
tered components are assumed to be equivalent in the 2P- 
HG function. For values of the relative amplitude parame- 
ter, c 2 , midway between totally forward or totally backward 
scattering, the properties of the forward and backward 
scattered components are dependent on each other. This 
is not the case with the 3P-HG function, where there are 
two half-width parameters, one for each scattering com- 
ponent. Until there are good disk-integrated reflectance 
measurements beyond 140° phase angle the nature of the 
scattering properties can only be estimated. Considering 
the phase angle coverage available in the phase curves 
presented here, the 2P-HG solutions are better con- 
strained. 

Anisotropic Multiple Scattering Results 

General. The RMS values of the anisotropic modeling 
results are within 0.5% of the isotropic modeling results, 
showing that statistically both fit the data equally well. 
However, the properties and shape of the single particle 
scattering function influence the values of the other param- 
eters in Hapke’s model. Below we discuss the different 
parameter values between the isotropic and anisotropic 
multiple scattering modeling results shown in Tables II and 
III in addition to examining the single particle scattering 
function predicted by both models. 


Single scattering albedo. The values of the single scat- 
tering albedo in the isotropic versus the anisotropic multi- 
ple scattering model results for Europa are within 5%. 
The isotropic model tends to produce slightly higher single 
scattering albedos. The modeling results for the darker 
surface of Ganymede show single scattering albedos within 
less than 1%, which are well within the error bars. These 
results demonstrate that for a surface as bright as Europa 
the two models predict qualitatively similar single scatter- 
ing albedos, but quantitatively the isotropic model pro- 
duces slightly brighter estimates. 

Roughness . Comparisons of the modeling results for 
the darker surface of Ganymede (as compared to Europa) 
show that both isotropic and anisotropic models predict 
the same macroscopic roughness properties to below 1% 
and well within the error bars. For the brighter surface of 
Europa, however, the differences in the predicted macro- 
scopic surface roughness are slightly larger. Both models 
predict a much smoother surface as compared to 
Ganymede or Callisto, and are thus qualitatively similar. 
The quantitative differences are with the limits of the 
parameter error bars (for the 0.55-^m Europa phase 
curve results the differences are at the edge of the error 
bars). 

Opposition effect. One of the largest differences ob- 
served between the isotropic and anisotropic multiple scat- 
tering models is between the opposition effect parameters 
(bo and h). Even so, the solutions to Ganymede’s phase 
curve have parameter values which are qualitatively similar 
and the differences are within the error bars for both oppo- 
sition parameters. This is also true for the solutions to both 
0.47- and 0.55-/xm leading phase curves of Europa. The 
phase curves for which the isotropic and anisotropic model- 
ing results show differences outside of the estimated error 
bars are the Europa trailing hemisphere phase curves. 

Single scattering function. Within the error bars the 
single particle scattering functions predicted by both iso- 
tropic and anisotropic models are the same for all analyses 
except for Europa’s leading hemisphere. The results are 
plotted in Figs. 11 and 12. Phase angle coverage for this 
hemisphere extends from essentially 0° to 120°. In this 
phase angle range the two models predict qualitatively 
similar single scattering behavior. The anisotropic solution 
predicts a slightly stronger backscattering component. At 
phase angles greater than 130° the two models drastically 
diverge. The isotropic multiple scattering model predicts 
a substantial forward scattering component, whereas the 
anisotropic multiple scattering model predicts nearly no 
forward scattering behavior. There are no data beyond 
120° phase angle that would constrain the models in the 
forward scattering regimes. 
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DISCUSSION 

Europa 

Europa has been observed to show a strong hemispheri- 
cal dichotomy in its rotational lightcurve (Stebbins 1927, 
Stebbins and Jacobsen 1928, Millis and Thompson 1975, 
Morrison and Morrison 1977, Johnson et al 1983, McEwen 
1986, Nelson et al 1986) where the leading hemisphere 
is brighter than the trailing hemisphere. McEwen (1986) 
showed that the hemispherical dichotomy could be mod- 
eled as a cosine function from the apex of orbital motion 
based on mosaics of the Voyager images, which is consis- 
tent with magnetospheric interactions with the surface 
(Pospieszalska and Johnson 1989). The hemispherical 
asymmetries on Europa have been attributed to two exo- 
genic processes: (1) preferential bombardment of the lead- 
ing hemisphere by micrometeorites (Shoemaker and Wolfe 
1982, Clark et al 1983, Clark and Lucey 1984) and (2) 
preferential sputtering by, and implantation of, magneto- 
spheric ions on the trailing hemisphere (Lane et al 1981, 
Clark et al 1983, Johnson et al 1984). 

The main spectral component of Europa’s surface is 
water ice (Pilcher et al 1972) with only a few percent 
impurities (Clark 1980). Clark (1980) modeled the spectral 
characteristics of the icy Galilean satellites and found that 
Europa’s spectrum most closely matches that of medium- 
to-fine grained water frost. Lane et al (1981) reported 
the detection of an absorption feature at 0.28 pm in the 
reflectance spectrum of Europa’s trailing hemisphere taken 
with IUE, which they identified as an S0 2 gas absorption 
feature. They hypothesized that this S0 2 feature was a 
product of sulfur ion implanatation into a water ice matrix, 
providing evidence for the magnetospheric alteration of 
Europa’s surface. More recently, however, Noll et al (1995) 
made spectroscopic observations of Europa using the Hub- 
ble space telescope and identified a broad absorption fea- 
ture at 0.28 p,m which they matched to S0 2 frost, which 
is consistent with the laboratory work of Sack et al (1992) 
and indicates a possible endogenic origin for the S0 2 seen 
on Europa’s trailing hemisphere. No clear detection of 
these absorption features has been made on Europa’s lead- 
ing hemisphere. Analysis of Europa’s UV and visible spec- 
tra indicates the presence of sulfur dioxide and sulfur com- 
pounds (Lane et al 1981, Sack et al. 1992, Calvin et al 
1995, Noll et al 1995, Spencer et al 1995); however, it is 
still debatable if the source is exogenic, endogenic, or a 
combination of these processes. Sack et al (1991, 1992) 
showed that the general reddening in the UV of Europa’s 
trailing hemisphere is not necessarily due to the implanta- 
tion of sulfur or the deposition of sulfur dioxide, but is 
caused by the morphological changes in the ice grain struc- 
ture from ion bombardment effects. 

Sack et al (1992) state that the ion bombardment of 
water ice causes an annealing effect which decreases the 


number of microcracks and bubbles which can act like 
internal scatterers. According to McGuire and Hapke 
(1995) and Hartman and Domingue (1997) particles with 
moderate to high densities of internal scatterers are more 
backscattering than those with little to no quantities of 
internal scatterers. Therefore single particle scattering 
functions which are more backscattering on the leading 
hemisphere than the trailing hemisphere are consistent 
with preferential ion bombardment on the trailing hemi- 
spheres. All modeling solutions presented in this paper 
show single particle scattering functions with backscatter- 
ing components that are larger for the leading hemisphere 
solutions than the trailing hemisphere solutions. The extent 
of the forward scattering component appears to vary with 
hemisphere and with model. The isotropic multiple scatter- 
ing 3P-HG model predicts a significant forward scattering 
component for the leading hemisphere and essentially little 
forward scattering on the trailing hemisphere. The aniso- 
tropic multiple scattering 3P-HG model, however, predicts 
essentially little forward scattering for both hemispheres. 
The isotropic and anisotropic 3P-HG models predict very 
similar scattering functions for Europa’s trailing hemi- 
sphere, whereas the isotropic 2P-HG model predicts a sig- 
nificant forward scattering component. The isotropic 2P- 
HG and anisotropic 3P-HG solutions are similar for the 
leading hemisphere scattering functions, whereas the iso- 
tropic 3P-HG model predicts significant forward scattering. 
However, for the range of phase angles for which data 
are available, all three models show the same scattering 
behavior to wihin the parameter error bars. For the phase 
angle range for which coverage is available all three model 
solutions are consistent with the spectral characteristics of 
Europa’s surface and the preferential ion bombardment 
of the trailing hemisphere. The deviations between the 
modeling results arise for those phase angles for which 
there are no data available to constrain the model parame- 
ters. Verbiscer and Veverka (1990), Verbiscer et al. (1990), 
and Domingue et al (1997) have studied the scattering 
functions of terrestrial snow and have shown that regard- 
less of the scattering function used in the modeling; terres- 
trial snows are forward scattering. Since water ice/frost is 
a main constituent of Europa’s regolith, it is reasonable 
to hypothesize that the average scattering function of Euro- 
pa’s regolith may contain a forward scattering component. 
However, until very large (>140°) phase angle data are 
available the forward scattering nature of Europa’s regolith 
may remain a mystery. 

The values of the single scattering albedo for all model 
solutions are consistent with the lightcurves; the leading 
hemisphere values in both wavelengths are brighter than 
the trailing hemisphere values. The single scattering albe- 
dos from the isotropic 2P-HG and isotropic 3P-HG models 
are the same within the error bars. The anisotropic 3P- 
HG model tends to predict slightly darker single scattering 
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albedos than the isotropic 3P-HG model; however, the 
values are qualitatively similar. The small differences 
(<0.03) are possibly a reflection of similarly small varia- 
tions in other parameter values. 

The low values of 0-bar (<15°) found by all models 
indicate that the small roughness seen on the 2-km scale 
by the Voyager images extends down to the centimeter 
and millimeter scales. To within the error bars, the isotropic 
2P-HG and 3P-HG models return similar values of the 
roughness parameter. The anisotropic 3P-HG model solu- 
tions show, in general, lower, and therefore smoother, sur- 
face roughness. One characteristic of the surface roughness 
parameter is that its value should be independent of wave- 
length; thus modeling solutions at 0.47 and 0.55 fx m should 
overlap within the roughness parameter error bar limits. 
This is true for the isotropic multiple scattering model 
solutions; however, this is not the case for the anisotropic 
3P-HG model solutions for Europa’s trailing hemisphere. 
The anisotropic modeling results for the surface roughness 
parameter solutions to the 0.47- and 0.55-fxm phase curves 
do not overlap within the error bar limits. 

The value of the opposition parameters indicates that 
both hemispheres are composed of a fairly transparent 
material (consistent with the 90% or greater water ice con- 
tent estimated by Clark (1980)); however, the isotropic 3P- 
HG solution values of bo indicate that the trailing hemi- 
sphere material is slightly more opaque, which is consistent 
with the observations of more dark material (possibly the 
sulfur and sulfur dioxide seen by Sack et ai (1992) and 
Noll et ai (1995)) being mixed with the water ice on the 
trailing hemisphere (sulfur dioxide frost is dark only in 
comparison to material with reflectance albedos of 80% 
or greater; Nelson et ai (1980)). The opposition width 
parameter from both isotropic multiple scattering models 
predicts a uniform particle size distribution and porosity 
on Europa’s surface, which is consistent with a balance of 
ion bombardment on the trailing hemisphere and microme- 
teorite bombardment on the leading hemisphere. This con- 
tradicts the results of Buratti et ai (1988) who found the 
trailing hemisphere to be more porous than the leading 
hemisphere based on the analysis of IUE observations. 
The difference between the Buratti et ai (1988) findings 
and those in this paper (and Domingue et ai (1991)) is 
due to the absence of opposition measurements in their 
UV phase curves. The anisotropic 3P-HG model results for 
the opposition parameters are consistent with the isotropic 
results except for the solution values for Europa’s 0.55-/um 
trailing phase curve. The anisotropic multiple scattering 
solutions for the opposition parameters for this phase curve 
are equivalent to a more transparent paticle in a more 
compact regolith. However, the opposition width parame- 
ter (h) measures surface properties which are independent 
of wavelength. Therefore model h values should be inde- 
pendent of wavelength. This is true for the isotropic 2P- 


HG and 3P-HG model solutions, but this is not the case 
for the anisotropic 3P-HG solutions for Europa’s trailing 
hemisphere. Therefore the remainder of this discussion will 
focus on the isotropic multiple scattering model solutions. 

Quantitative estimates of the surface porosity on Europa 
can be made assuming a particle grain size distribution 
function. The opposition width parameter h is related to 
porosity by 

/! = -(») In (P)y, 

where P is the porosity (percentage of void space per unit 
volume) and Y is the grain size distribution function. If a 
grain size distribution similar to the Moon’s is assumed, 
then y is given by 

Y = V3 / In (r,/r s ), 

where r { is the radius of the largest particle and r s is the 
radius of the smallest particle in the size distribution 
(Hapke 1986). Figure 14 plots the porosity of Europa’s 
surface as a function of rj/r s . The porosities of the other 
icy Galilean satellites made using the same assumptions 
are plotted for comparison along with the results from 
Hartman and Domingue (1997) for the Moon. The high 
porosity for Europa is consistent with a medium to fine- 
grained water frost. 

Ganymede 

Ganymede also displays a hemispherical dichotomy in 
its visible lightcurve which is in the same sense as Europa 
but not as pronounced (Stebbins 1927, Stebbins and Jacob- 
sen 1928, Millis and Thompson 1975, Morrison and Mor- 
rison 1977, Johnson et ai 1983). The exogenic processes 
at work on Europa are also proposed to be the sources of 
the hemispherical differences on Ganymede. 

Water ice has also been identified on the surface of 
Ganymede (Pilcher et ai 1972), with a larger fraction of 
impurities present than estimated on Europa. If the water 
ice is intimately mixed with the darker non-ice component, 
then the estimated water ice abundance is on the order of 
90 wt% (Clark 1980); however, areal mixing models show 
that the water ice coverage on Ganymede could be as low 
as 50% and still agree with the spectral data (Spencer 1987). 
The best fits to Ganymede’s spectra show that the water 
absorption bands are best modeled as frost on ice (Clark 
1980) with the frost grain size larger than those predicted 
for Europa (Clark 1980, Clark et ai 1986, Calvin et ai 
1995). The spectra also indicate that the grain sizes on 
Ganymede's trailing hemisphere are larger than those on 
the leading hemisphere (Clark et ai 1986, Calvin et ai 
1995). In comparing Ganymede’s spectra with Europa’s, 
Nelson et ai (1987) noticed an additional source of absorp- 
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FIG. 14. Porosity profiles for Europa (dashed line), Ganymede’s leading (triangles) and trailing (squares) hemispheres, Callisto's leading (circles) 
and trailing (crosses) hemispheres, and the lunar near side (diamonds). These porosity values are given as a function of the ratio of the radii of 
the largest to smallest particle of the regolith assuming a grain size distribution typical of comminution processes. 


tion in Ganymede’s trailing hemisphere spectra shortward 
of 0.28 /xm which they attributed to the presence of molecu- 
lar oxygen in the form of ozone. Noll et ai ( 1 996) observed 
both hemispheres of Ganymede in the U V with the Hubble 
space telescope and detected the presence of an absorption 
band on Ganymede’s trailing hemisphere at 0.26 /xm which 
they identified as the Hartley band of ozone. This absorp- 
tion band was not present on the leading hemisphere. Both 
Nelson et al. (1987) and Noll et ai (1996) attribute the 
presence of ozone on the trailing hemisphere to ion bom- 
bardment and implantation process. Spencer et ai (1995) 
and Calvin et ai (1995) report the detection of absorption 
features in the visible at 0.577 and 0.627 /xm due to the 
presence of solid diatomic oxygen on Ganymede’s trailing 
hemisphere, which is consistent with the preferential inter- 
action of the magnetosphere with the trailing hemisphere. 

All model solutions are consistent with the spectral ob- 
servations of Ganymede’s surface structure and character- 
istics. They all show that the single scattering albedo is 
darker on Ganymede than Europa, which is consistent with 
geometric albedo measurements of these two satellites. 
Both isotropic multiple scattering models find that the sin- 
gle scattering albedo at 0.55 /xm is brighter on the leading 
hemisphere than the trailing hemisphere (which is consis- 


tent with the visible lightcurves); however both models 
show that the opposite is true at 0.47 /xm, which is unex- 
pected. The anisotropic multiple scattering model was not 
fit to the 0.55-/xm data (due to lack of phase angle data at 
larger, >90°, phase angles); therefore similar comparisons 
cannot be made. A wavelength dependence is seen in some 
of the model solutions. The trailing hemisphere at 0.47 /xm 
has a nearly isotropic scattering function predicted by the 
models, whereas at 0.55 /xm the solutions show a strongly 
backscattering behavior. Each model predicts a backscat- 
tering behavior for Ganymede’s leading hemisphere at 0.47 
/xm. However, the isotropic 2P-HG model predicts a for- 
ward scattered component at 0.55 /xm. 

The roughness dichotomy measured by all models agrees 
within their error bars. A rougher trailing hemisphere is 
also consistent with a lower physical albedo on the trailing 
hemisphere. It is noted that the 0.55-/xm phase curves for 
both hemispheres of Ganymede lack the necessary large 
phase angle observations to make realistic determinations 
of the Hapke roughness parameter. However, the agree- 
ment with the 0.47-/xm Hapke roughness parameters lends 
validity to the 0.55-/xm roughness measurements since the 
roughness parameter does not vary with wavelength. 

The opposition parameter values of all the models are 
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consistent with the evidence that there is a larger percent- 
age of non-ice component on Ganymede compared to Eu- 
ropa. The bo values for Ganymede at 0.47 ptm indicate 
that more opaque material is present on Ganymede than 
on Europa; however, the Q.55-(jLm results give conflicting 
measurements. The spectral work of Clark (1980), Clark 
et al (1986), and Calvin et al (1955) estimates that the 
particle grain sizes on Ganymede are generally larger than 
those on Europa, smaller than those on Callisto, and have 
the added complication that they are larger on the trailing 
hemisphere than the leading hemisphere. The opposition 
width parameter is physically related to both porosity and 
grain size distribution. The smaller values of h are generally 
correlated to larger porosities or small compaction states; 
however the correlation with grain size is less straightfor- 
ward since it is more strongly dependent on the grain size 
distribution (the ratio of largest grain to smallest grain). 
If smaller values of h are correlated with smaller grain 
sizes, then the spectral analyses would predict that the 
values of h in increasing order would be Europa's surface 
(and possibly the leading hemisphere of Callisto), the lead- 
ing hemisphere of Ganymede, the trailing hemisphere of 
Ganymede, and the trailing hemisphere of Callisto. The 
Hapke model solutions show that the only h values which 
do not follow this trend are the values for Ganymede's 
trailing hemisphere and Callisto's leading hemisphere. This 
trend, however, assumes similar porosities and that smaller 
grain sizes correspond to smaller grain size distributions, 
which are extremely unlikely. The plots in Fig. 14 compare 
the porosity estimates from the isotropic multiple scatter- 
ing solutions of both hemispheres of Ganymede assuming 
a lunar-like particle distribution. The anisotropic multiple 
scattering solutions are similar to the isotropic multiple 
scattering modeling results and are, therefore, not in- 
cluded. The porosity of the leading hemisphere is very 
similar to Europa’s surface; however, the trailing hemi- 
sphere is more compact than the lunar regolith. This com- 
paction state is highly unlikely in view of the match be- 
tween Ganymede's spectra and that of frost on ice. This 
indicates that the grain size distribution function assumed 
in the derivation of these porosities is not indicative of the 
grain size distribution on the surface of Ganymede. 

Callisto 

Callisto displays a hemispherical dichotomy in its visible 
lightcurve which is in the opposite sense to that seen on 
Europa and Ganymede (Stebbins 1927, Stebbins and Ja- 
cobsen 1928, Millis and Thompson 1975, Morrison and 
Morrison 1977, Johnson et al. 1983); i.e., the trailing hemi- 
sphere is brighter than the leading hemisphere. Buratti 
(1991) and Calvin and Clark (1993) conclude from their 
analyses that the dominant process affecting the optical 


surface of Callisto is micrometeoritic gardening. Meteorite 
bombardment is thought to be the dominant exogenic pro- 
cess at work on Callisto since its orbit is distant enough 
to minimize magnetospheric interaction effects. 

Water ice absorptions have also been detected on Cal- 
listo (Pilcher et al. 1972). However, it has the largest 
amount of dark non-ice component incorporated into its 
surface among the three icy Galilean satellites. Clark 
(1980) estimated that the amount of water on Callisto 
ranges from 30 to 90 wt% assuming an intimate mixture 
of ice and non-ice component. In contrast, modeling studies 
by Calvin and Clark (1991) and Roush et al (1990) showed 
that the weight percent of water on Callisto could be as 
low as 20 to 45% if combinations of intimate and areal 
mixtures of the water ice and non-ice components are 
considered. Spencer (1987) showed that areal mixtures of 
10% water ice are consistent with the spectral data. The 
water ice grain size on Callisto is the largest of the three 
satellites, being predominately coarse grained (Calvin and 
Clark 1991, Roush et al 1990). Calvin and Clark (1993) 
and Calvin et al (1995) noted, however, that there is a 
water ice grain size dichotomy on Callisto, with the spectral 
signature of fine grained water ice in the 3- to 3.6-/um 
region matching that of the leading hemisphere. 

The largest differences in surface structure predicted 
by the isotropic 2P-HG and 3P-HG models are for the 
surface of Callisto. The values of the single scattering 
albedo found with the isotropic 3P-HG model are consis- 
tent with the visible lightcurves; the trailing hemisphere 
values in both wavelengths are higher than the leading 
hemisphere values. The isotropic 2P-HG model shows 
the opposite trend which is inconsistent with the 
lightcurves. If the observations of a more dark, non-ice 
component on Callisto than on the other two icy Galilean 
satellites are accurate, then the bo values returned by 
the models should be higher for Callisto, indicative of 
the presence of more opaque material. Both models are 
consistent with this predicted trend with the anomalous 
exception of the 0.47-/xm solutions for trailing hemi- 
sphere. The values of the opposition width parameter, 
h , are consistent with a more porous trailing hemisphere. 
The spectral evidence suggests that the grain sizes on 
the trailing hemisphere are larger than those on the 
leading hemisphere, which may be similar in grain size 
to that seen in Europa’s surface. This would imply a 
grain size distribution on the trailing hemisphere smaller 
than that on the leading hemisphere, thus raising the 
possibility of a higher small grain frost content on the 
leading hemisphere. The porosities plotted in Fig. 14 are 
consistent with this scenario. Without any large phase 
angle data for the leading hemisphere any discussion of 
hemispherical differences in surface roughness or single 
particle scattering function is speculation only. 
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CONCLUSIONS 

The contrasts between our results and those of Buratti 
(1991) for Ganymede and Callisto indicate that there is a 
vast difference in the predicted physical characteristics of 
a planetary regolith dependent on the choice of single 
particle scattering function incorporated in Hapke's model. 
The results presented here show that the general physical 
characteristics predicted by a Hapke model (with isotropic 
multiple scattering) incorporating a 2P-HG scattering func- 
tion are similar to the predictions given by a similar Hapke 
model incorporating a 3P-HG function as long as the data 
set being modeled has adequate coverage in phase angle. 
Conflicting results occur when the large angle coverage 
is inadequate. For surfaces as bright or darker than 
Ganymede both isotropic and anisotropic models predict 
similar photometric behavior over the phase angle ranges 
of the data set. For phase curves with phase angle coverage 
less than 120° the isotropic and anisotropic models predict 
similar results (within that limited phase angle range) for 
objects as bright as Europa; however, they may predict 
different single scattering behavior for phase angle ranges 
not included in or constrained by observations. The photo- 
metric surface textures predicted in our results support the 
spectral evidence for magnetospheric alterations of the 
optical properties of the surface of Europa and Ganymede. 
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